Destruction of chain-superconductivity in YBa 2 Cu 4 08 in a weak magnetic field 
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We report measurements of the temperature dependent components of the magnetic penetration 
depth A(T) in single crystal samples of YBa2Cu40s using a radio frequency tunnel diode oscillator 
technique. We observe a downturn in A(T) at low temperatures for currents flowing along the b and 
c axes but not along the a axis. The downturn in Xt is suppressed by a small dc field of ~ 0.25 T. 
This and the zero field anisotropy of A(T) likely result from proximity induced superconducting on 
the CuO chains, however we also discuss the possibility that a significant part of the anisotropy 
might originate from the CuC>2 planes. 



The Y-based high-T c cuprate superconductors are 
unique in that they have quasi-one-dimensional CuO 
chain structures in addition to the Cu02 planes in which 
the interactions responsible for superconductivity are 
thought to originate. The effect of these chain layers 
on the normal and superconducting state properties has 
long been debated. Indeed, since the Y-based cuprates 
are so far the only hole-doped cuprates which show quan- 
tum oscillations in their underdoped state,— ~— it is natural 
to question if the coherent electron orbits could in fact 
originate in the chain parts of the Fermi surface (FS).— £ 
A important question is whether the chain FS is super- 
conducting and if so is this superconductivity quenched 
by a much smaller magnetic field than the plane. 

YBa2Cu307_,5 (Y-123) contains a single CuO chain 
per unit cell whose oxygen content can be varied 
to tune the doping level on the Cu02 plane. Its 
close relative YBa2Cu40s (Y-124) has a double chain 
layer that is stoichiometric and a planar state that 
is slightly underdoped.— According to band-structure 
calculations^^ the chains form quasi- ID electronic 
bands which cut across the Cu02 plane FS sheets, with 
considerable hybridization between the states close to the 
crossing points. As the electron-pairing interaction which 
gives rise to superconductivity in the cuprates is likely 
to rely strongly on the electronic structure of the quasi- 
two-dimensional plane, it might seem unlikely that there 
any intrinsic pairing of electrons on the chain. Indeed, in 
isostructural Pr-124, the double chain network is metallic 
yet exhibits no superconductivity down to 0.5 

In both Y-123 and Y-124, the CuO chains have been 
modelled as an essential normal layer which is coupled 
to the planes via single electron tunnelling, similar to 
the classical proximity effect between normal metals and 
superconductors This is augmented by the strong hy- 
bridization between plane and chain states which occurs 
at certain momentum values. This model predicts that 
the chain states will have low energy gap structures aris- 
ing from variations of the gap within the chain FS. This 
implies that the superfluid density should have quite dif- 
ferent temperature dependencies for screening currents 
flowing in the a or b (chain) directions. Experimentally 



however, it is found that for Y-123 the temperature de- 
pendence of A and Xb are very similar 11 ' 12 although their 
zero temperature values differi 13 ' 14 This has led to the 
proposal that there are intrinsic pairing interactions on 
the chains and that the planes and chains are predom- 
inately coupled by Josephson-like pair tunnelling*^ Al- 
ternatively, this has been explained by chain disorder. 15 

In this paper, we report measurements of the 
anisotropic components of the magnetic penetration 
depth in Y-124 single crystals as a function of tempera- 
ture and dc field. A rapid suppression of the supercon- 
ducting component on the CuO double chain is observed 
in a small applied field. These results are significantly dif- 
ferent to those obtained on optimally doped Y-123 and 
appear to support the proximity effect model for Y-124. 
However, some questions remain. 

Single crystals of YP^CiuOs were grown using a high 
oxygen pressure flux based method^ Penetration depth 
was measured using a radio frequency (RF) tunnel diode 
oscillator method^ operating at ~ 12 MHz. The sample 
is attached with vacuum grease to a sapphire rod and is 
placed in a copper coil which forms part of the oscillator's 
tank circuit. A small superconducting solenoid allows us 
to apply a dc field co-linear with the weak RF probe field 
(7J RF ~10- 6 T). 

Changes in the resonant frequency of our RF oscil- 
lator AF are directly proportional to changes in sam- 
ple's superconducting volume as a function of temper- 
ature or field. The samples are thin platelets with di- 
mensions £ a ,b,c along the respective crystallographic di- 
rections (Fig. Q]). With the RF field applied along the b 
direction the change in frequency AF due to changes in 
the A Q and A c is given by AF = 2a/3(£ a £ b AX a + £ b £ c AX c ) 
where a is a constant set by the coil geometry and /3 is 
the effective sample demagnetizing factor (here we ne- 
glect the small contribution from the currents on small 
sample surfaces perpendicular to the field). Hence, con- 
tributions from two components of A are always mixed. 
In principle, it is possible to separate the components 
by making measurements of A A with TJrf in ah three 
directions, however in practice this approach is inaccu- 
rate because of uncertainties in the demagnetizing fac- 
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FIG. 1: (color online) Temperature dependence of the penetration depth for the three principle directions for YBa2Cu4C>8. 
The shapes of the samples are shown in the Aj, panel. Their dimensions (£ a X lb x £ c ) are: sample 1: 115 x 80 x 15/xm 3 , sample 
2: 60 x 420 x 8 pm 3 , sample 3: 170 x 60 x 6 pm 3 . The RF susceptibility (x') close to T c (normalized to —1 at low temperature) 
is shown in the A c panel. 



tor particularly in the H\\c configuration. Instead, we 
cleaved the sample (using a razor blade) halving l a and 
thus doubling the contribution of AA C J^ Measurements 
of sample 1 with H\b (where the demagnetising effect is 
small) before and after cleaving thus allow us to isolate 
the A c contribution which then can be subtracted from 
the measurements with H \\a and H \\b to yield the three 
components of A. Consistency was checked by cleaving 
the sample a second time. The extracted AA C (T) was 
used to extract AA a and AA& from two further samples. 

The temperature dependence of the three extracted 
components of A below T = 40 K are shown in Fig.[TJ The 
data for AA a (T) and AAf,(T) for all three samples is very 
consistent. A a has a linear T dependence from base tem- 
perature up to ~ 15K, with slope d\ a /dT = 9.6(6)A/K 
which is consistent with a simple ci-wave model. The 
behavior of AAf, and AA C are quite different to AA a . 
Between ~ 15 and 30 K AA&(T) is linear with a slope 
dXb/dT = 24(2)A/K which is ~2.5 times larger than for 
A a . Below ~15 K however, there is a downturn in AAb(T) 
which indicates the onset of additional screening which 
reduces A& by ~150A at T = 2K compared to the ex- 
trapolated linear behavior. Similar behavior is found for 
A A c where the downturn sets in at approximately the 
same temperature and d\ c /dT = 130A/K. A downturn 
in X(T) had been seen in previous measurements of Y- 
12 4 18 i 19 however the separate contributions of A a and Xb 
were not determined. 



In the conventional model of <i-wave 
superconductivity 2 ^ with a circular 2D FS we ex- 
pect = ^A/i^i^d • W e cannot measure A(0) directly 
in our experiments. However, infra-red reflectivity 
(IRR) experiments^ on Y-124 give A a (0) = 2000 A and 
A 6 (0) = 800 A, (the anisotropy A a /A fc = 2.5). Al- 
though the A anisotropy measured by IRR for Y-123 
{X a /Xb — 1.6) is somewhat higher than that measured 
by other other techniques (A a /Ab — 1.2)^ it does seem 
likely that X a /Xb is signihcantly greater than unity for 
Y-124 which is in the opposite direction to our measured 
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FIG. 2: (color online) Normalized superfluid density [p — 
(A(0) 2 /(A(0) + AA(T)) 2 ] for sample 2 in all three directions 
using values of A(0) taken from infra-red measurements^ 3 " 
(A a (0) = 2000 A and A 6 (0) = 800A) and aligned polycrys- 
talline measurements (A c (0) = 6150 A).— The inset shows 
theoretical predictions for p in the proximity coupling (single 
electron tunneling) model taken from RefUCI. 



anisotropy in dX/dT. Hence, Fermi velocity anisotropy 
alone cannot explain the data within this simple model. 
If the order parameter had a significant s component 
(i.e., d + rjs) so that the nodes move away from <f> = 45° 
towards the 6-direction, then the paramagnetic current 
produced by the thermally excited nodal quasiparticlcs 
could be larger in the 6-direction and this, in principle, 
could overcome the opposite anisotropy of A(0). Within 
a single band model a value of X a /Xb — 2.5 implies a 
substantial anisotropy of the in-plane Fermi velocity. 
Calculations based on fits to photoemission results 21 
however, suggest that at a doping level of p = 0.14 the 
anisotropy of the in-plane superfluid density is only 
- 2%. 

An explanation for both the larger value of A(0) and 
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AA(T) in the 6-direction can be found in the plane-chain 
proximity models mentioned above. In Fig. [5] we show the 
normalized superfluid density p = A(0) 2 /(A(0) + AA(T)) 2 
for sample 2, calculated using the values of A(0) taken 
from infra-red measurements^ 3 - (A a and Afc) and aligned 
polycrystalline measurements (A c ).— The behavior in the 
a-direction is similar to that found for p a and pb of Y- 
123, and varies linearly with T up to ~ T c /2. How- 
ever, both pb{T) and p c (T) are quite different to p a (T) 
with both showing strong upward curvature for temper- 
ature below ~ T c /3. Note that the small downturn in 
AA(,(T) below 15 K noted above (Fig. [TJ) is not the pri- 
mary reason for the disparity between p a and pb in Y- 
124; rather it is caused by the substantial difference in 
the ratio (dX/dT)/X(0) which is 4.8 x lO^K" 1 for the a 
axis but ~ 6 times larger (30 x 10~ 3 K^ 1 ) for the b axis. 
Clearly small uncertainties in the values of A(0) will not 
have a large effect on this. These upturns are similar to 
those reported for the average in-plane and c axis super- 
fluid density measured on polycrystalline samples* 2 ^ 

In the inset to Fig. [2] we show the predictions of the 
single electron tunnelling plane-chain proximity coupling 
model of Reffiol (similar results using a similar model 
are also given in RefQ). These predictions are very sim- 
ilar to our experimental findings, with both Pb{T) and 
p c (T) having strong upward curvature below ~ T c /3. In 
this particular simulation both plane and chain FS were 
assumed to be superconducting, but similar results are 
found when there is no intrinsic gap on the chains.— A 
feature of this proximity model is that the upward cur- 
vature of pb is strongly suppressed by impurity scattering 
(e.g. broken chain segments)^ The difference between 
Y-123 and Y-124 may then be due to fact that Y-124 
has completely full chains whereas in Y-123 the chain 
segments are much shorter (more disorder). Inter-chain 
coupling in Y-124 will reduce localization effects which 
are likely to be present in the more ID Y-123 chains. 

Next we discuss the field dependence of A. Here we 
apply a dc field co-linear with the RF probe field along 
the c direction where H c i is maximal. Data for sam- 
ple 2 are shown in Fig. [3] for dc fields up to 50 mT. In 
this field configuration the measured AA a & is a mixture 
of AA a and AAf, but because of the aspect ratio of this 
particular sample {lbji a — 5) around 80% of AA & comes 
from AXb- In this configuration the calibration factor is 
difficult to calculate (especially for a sample with such 
an elongated shape), so we determined it experimentally 
by measuring, in both the H\\ab and H\\c configurations, 
a twinned sample of Y-123 with almost identical shape. 
From the total frequency shifts measured when the sam- 
ple was withdrawn from the coil and from the sample 
dimensions we estimate that the demagnetizing field en- 
hancement factor j3 ~ 5, hence the effective surface fields 
are 5 times larger than the applied field. 

For ffdc = (Fig. [T]) a strong downturn in X a b{T) 
is seen below T ~ 15K as in the direct measurements of 
AXb(T) (Fig.[T]). As i?dc is increased, initially there is no 
change and then at poH^c — 5mT, Xb(T = 2K) starts 



to increase monotonically and eventually saturates for 
PoHdc ^ 50 mT (see inset Fig. |3|). The temperature de- 
pendence of X a b also changes, from ~ T 7 for H^ c = 
to ~ T 1 - 2 for /io-ffdc = 50 mT. The applied dc field wipes 
out the downturn in A a b and returns the temperature 
dependence to a quasi-linear variation expected from a 
conventional d-wave superconductor with small amounts 
of impurities^ Similar behavior with comparable field 
scales was observed in two other samples. The high field 
behavior of X a b(T) is less linear than A a at the lowest tem- 
peratures perhaps because of anisotropic impurity scat- 
tering. 

For poHdc ^ 35 mT AA & is perfectly reversible on 
cycling from base temperature (where the field was in- 
creased) up to 20 K and back again. For fields above 
this AA a b(T) was hysteretic on the first temperature cy- 
cle and then reversible thereafter. This irreversibility in- 
creased linearly with field (see inset Fig. [3]). We attribute 
this hysteresis to vortex motion. Initially, as the vortices 
enter the sample they are not in deep pinning wells and 
so give a strong contribution to the measured effective A. 
Upon temperature cycling they settle on strong pinning 
sites (or move to the center of the sample) and therefore 
no longer give a strong contribution to A. As the field at 
which X a b{T) begins to change is much below the field at 
which the hysteretic behavior set in, and since for fields 
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FIG. 3: (color online) Field dependence of X a t for sample 2 
(dimensions 60 x 300 x 8pm 3 ) from to 50 mT in 5mT incre- 
ments. The b-axis dimension is less than that stated in Fig. 
[T] because it was broken between runs. The aspect ratio of 
this sample means that ~80% of \ a b comes from Xt- To re- 
move artifacts arising from the field dependent (but tempera- 
ture independent) background from the measurement coil the 
data have been shifted in frequency so that they coincide at 
T = 20 K for all fields. The lower inset shows the change in 
A at the fixed temperature of 2 K along with the irreversible 
changes to A, showing the onset of flux entry at ~ 35 mT. 
The upper inset shows the superfluid density along a and b 
axes in zero field together with the b axis data in 50 mT. Here 
AAfc(T) at 50 mT has been approximated by 1.2AA a f, . The 
factor 1.2 accounts for the small contribution of A a to \ ab . 
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greater than this the reversible part of AA(if) begins to 
saturate whereas the irreversible part increases linearly 
with H , it is clear that vortex motion is not responsible 
for the changes in the reversible part of AX a b(H) which 
are shown in the main part of Fig. [3] Rather reversible 
changes reflect a decrease in the Meissner screening cur- 
rent. Direct measurements of the dc magnetization using 
a SQUID magnetometer showed that the field of first flux 
penetration was ^35 mT, for this sample and field con- 
figuration at T = 10 K, which coincides with the onset of 
irreversibility. 

Field dependence of A in the Meissner state can re- 
sult from several effects. A linear increase of A with 
H is expected from Doppler shifted quasiparticles close 
to the nodal points.— This effect is however, very small 
(~2A for fioH — 10 mT for Y-123) and has never been 
conclusively observed (the measured X(H) in Y-123 does 
not have the temperature or field orientation dependence 
expected from theory and may therefore have another 
origin) j 12 i 24 A decrease in A with increasing field can also 
result from the presence of Andreev bound states^ but 
this is in the opposite direction to what is observed here. 

Instead, the most likely explanation is that in Y-124 
the small dc field quenches a least part of the supercon- 
ducting current in the chains. We might expect this to 
occur when the Doppler shift of the chain quasiparticle 
energies are of order the chain gap. It is often observed 
that proximity-driven superconductivity is quenched in 
small fields^ for example the ir band superconductivity 
in MgB2^ In Y-124 the magnetic field does not affect 
A(T) above ~ 15 K and the strong anisotropy in dX/dT 
between the a and b direction as discussed above remains. 
Also the field induced change in X a b ~ 190 A at the low- 
est temperature which does not significantly change the 
anisotropy of A. This is illustrated in the inset to Fig. 
[3] where we show the field suppressed superfluid density 



along the b axis along with the zero field results. 

The fact that only part of the excess b axis super- 
fluid (relative to the a axis) is suppressed in weak field 
(^0.25 T including demagnetizing effects) might suggest 
that there is significant variation of the gap along the 
chain FS. 28 If the system is sufficiently clean, so that 
scattering between regions of the chain FS with differ- 
ent gaps (resulting from a variation of the plane-chain 
hybridization) is small, then the downturn in A&(T) and 
A C (T) for T <15K could result from a region of chain 
FS with weak pairing becoming superconducting. This 
region then becomes normal when a weak field is ap- 
plied. Alternatively, the downturn and field depen- 
dence could reflect the superconductivity in the chain 
as whole, and then the remaining anisotropy in high field 
would reflect the anisotropy of the CuC>2 planes. Al- 
though, as mentioned above, there is no indication of such 
anisotropy from photoemission experiments we note that 
quantum oscillation 2 ' 5 and magnetotransport2£ measure- 
ments show that there is a significant reconstruction of 
the FS in this material at low temperature. Recently, 
evidence of significant electronic anisotropy in the CuC>2 
planes of underdoped Y-123 has been found in Nernst ef- 
fect measurements^ and it is likely that properties of 
Y-124 are similar. Further theoretical work calculat- 
ing the effect of magnetic field on the superfluid den- 
sity in the proximity model a 9 ' 10 would help decide be- 
tween these two competing interpretations. Completing 
the picture of plane-chain superconductivity in Y-124 and 
Y-123 should be help our understanding of quantum os- 
cillations and the underlying electronic structure of these 
materials. 
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